Quantum and classical diffusion in small-world networks 
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We study numerically quantum diffusion of a particle on small-world networks by integrating the time- 
dependent Schrodinger equation with a localized initial state. The participation ratio, which corresponds to the 
number of visited sites in the case of classical diffusion, as a function of time is measured and the corresponding 
diffusion time r is computed. In a local regular network, i.e., in the network with the rewiring probability p — 0, 
the diffusion time depends on the network size N as r ~ N, while the behavior r ~ log N is observed as p 
becomes finite. Such fast diffusion of a particle on a complex network suggests that the small-world transition 
is also the fast-world transition from a dynamic point of view. The classical diffusion behavior is also studied 
and compared with the quantum behavior. 
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There has been a surge of research activity on various as- 
pects of complex networks since some important features of 
real networks were successfully explained by simple model 
networksii 2 * 3 - In particular, the Watts -Strogatz (WS) model 2 
was the first to produce networks with small-world behavior, 
characterized by the path length increasing logarithmically 
with the network size. Subsequently, a number of studies have 
been performed on such complex networks, focused mostly on 
structural properties of the networks. On the other hand, ver- 
tices of a real network may have some internal degrees of free- 
dom, interwoven with the structure of the network. Motivated 
by this observation, a group of papers have also investigated 
statistical mechanical models defined on complex networks. 
For example, spin models like the Ising model and the XY 
model on WS networks have been shown to undergo finite- 
temperature phase transitions of the mean-field nature^ dis- 
closing the role of long-range interactions along the shortcuts. 
In some cases the time-development of the network structure 
may be coupled to the dynamics of degrees of freedom defined 
on vertices. Still the study of vertex dynamics without consid- 
ering its influence on the network structure can be useful as 
a first step toward the complete understanding. In this spirit, 
dynamic models defined on networks draw much attention: 
Epidemic spreading as well as classical diffusion on complex 
networks have been studied, 6,7 and very recently, the dynamic 
universality class of the XY model on the WS network has 
been identified. 8 

Properties of a quantum mechanical model put on a com- 
plex network have also been studied with respect to the spec- 
tral properties of the Laplacian operator 9 and the localization- 
delocalization transition in the presence of site disorder^ In 
this work we consider quantum as well as classical diffusion 
of a particle on the WS network without site disorder, and in- 
vestigate how the diffusion time scales with the network size. 
It is well known that as soon as the WS network (of size 7Y) 
has a finite fraction of shortcuts, it undergoes the small-world 
transition that the characteristic path length I changes its be- 
havior from I ~ N to log iVi 2 *!! ^From the dynamic point of 
view, this change from the large-world to the small-world be- 
havior is expected to be accompanied by a sharp change in the 
behavior of the diffusion time: If the world is small, the trav- 



eling time around the world should also be short. Indeed the 
diffusion time r associated with the participation ratio is ob- 
served to change the size dependence from r ~ N to log N, 
which is to be compared with the case of classical diffusion: 
r ~ ./V 2 to N. 

The low-dimensional tight-binding electron system has 
been studied extensively in relation to the localization transi- 
tion in the presence of disorder. In the absence of disorder, the 
energy eigenstate of a tight-binding electron on a local regular 
lattice is always in the extended state due to the Bloch theo- 
rem. Here we consider the tight-binding electron on the WS 
network, described by the time-dependent Schrodinger equa- 
tion 



ff|*), 



(1) 



where we have set H = 1 and the ket |^) has the position rep- 
resentation ^„ = (n|\&) at the nth vertex. The Hamiltonian 
in the position representation takes the form 



Hn 



A for n' e A,, 
otherwise, 



(2) 



where the on-site energy has been assumed to be uniform and 
set equal to zero (H nn = 0), A is the hopping energy, and A„ 
represents the set of neighbors of vertex n. For example, in 
a local regular network with the connection range r = 1, we 
have A„ = {n— 1, n+1}. 

The WS network is constructed according to the standard 
procedure in Ref. 2: A one-dimensional (ID) local regular 
network with the connection range r is built first, then with the 
rewiring probability p each local edge is rewired to a randomly 
chosen other vertex. (See Fig.^ illustrating the case r = 2.) 
Once the WS network is constructed in this way, we normalize 
the time t in units of 1/ A, and integrate the time-dependent 
Schrodinger equation given by Eqs. Q and l|2) numerically 
by means of the fourth order Runge-Kutta method with the 
discrete time step St = 0.01, starting from the initial condition 
that the electron is localized at randomly chosen vertex m, i.e., 
^„(t=0) = S n>m , For simplicity, we set the initial position 
m = 0. 
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FIG. 1 : Construction of the Watts-Strogatz (WS) network. Starting 
from the regular one-dimensional lattice in (a) with the connection 
range r — 2, we visit each edge and rewire it with given probability 
p, thus generating shortcuts in (b). Via this procedure a complicated 
network structure emerges, as shown in (c). 



In the absence of shortcuts, the analytic solution of Eq. Q 
is easily found: 12 In the simplest case of r = 1, the Fourier 
transformation = J2 n 1 &ne mk yields 



(e lfc +e- tk )^ k = 2cosfc* fe , 



which in turn leads to 

*fc(t) = * fc (0)e- 24 * cosfc . 
The inverse Fourier transformation 



(3) 



(4) 



(5) 



in the limit of N — > oo, combined with the Jacobi- Anger ex- 
pansion e lzcos9 = J2m=-N* m J m ( z ) eZme gi ves me solution 



|tf„(t)| = \J n (2t)\ 



(6) 



with the Bessel function J n of the order n. For 
r = 2, it is straightforward to obtain *f? n (t) = 

J2 p (-l) n - p Jn-2 P (2t)J p (2t). 

As time proceeds, the initial localized state diffuses and 
eventually evolves to an extended state in which ^ 
at any n. Figure [2] shows such diffusion of the wave packet 
on the WS network with the connection range r = 2 for the 
rewiring probability p = and 0.2. It is clearly shown that 
diffusion occurs much faster in (b), namely, for p ^ 0. In this 
case as soon as the wave packet diffuses to a nearby end point 
of a shortcut, the other end of the shortcut becomes a new cen- 
ter for diffusion at the next time step, resulting in faster spread 
of the wave packet. 

To describe the diffusion of the wave packet, one may 
measure the variance of the wave packet, 13 cr 2 (t) = 
{Sb(t)\x 2 \i&(t)} - (t)\x\^ (t)) 2 with the one-dimensional 
position x on the network. However, the variance a 2 tends to 
overestimate the diffusion of the wave packet. As an exam- 
ple, consider a state where the wave packet is localized at two 
vertices separated by a distance O(N). Although the state is 
well localized in the sense that the particle can be detected 
only at a few number of sites, the variance has a very large 
value 0(N 2 ). This is in contrast with a local network, where 
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FIG. 2: Spread of the wave packet initially localized at x — on 
the Watts-Strogatz network with the connection range r = 2 for the 
rewiring probability p = (a) and (b) 0.2. In both (a) and (b), the 
three curves described the wave packet at time t = 1.0, 5.0, and 10.0, 
respectively (from top to bottom). It is observed that the wave packet 
diffuses much faster in (b) than in (a). For clarity, 2|^| 2 (instead of 
|*| 2 for t = 1.0) is plotted for t = 5.0 and 10.0. 



the particle can hop to only locally connected nearby sites and 
thus a 2 may be used as a measure of the diffusion without 
any confusion. To avoid this difficulty with a 2 , we instead 
consider the participation ratio Pq (t) as a function of time: 



p (f) = En\*n(t)\ 2 



(7) 



which has the value unity for a wave packet localized com- 
pletely at one site. For a completely extended state, on 
the other hand, we have l^n] 2 ~ l/N and accordingly, 
the participation ratio Pq ~ N. The initial localized state 
with Pgit—O) = 1 thus evolves to the extended state with 
Pq(<— voo) = O(N) as the time-dependent Schrodinger 
equation is integrated in time. 

With Eq. (|6} and the asymptotic form of the Bessel func- 
tion, a lengthy but straightforward calculation leads to the 
asymptotic behavior Pq(t) ~ t on the ID regular net- 
work (without shortcuts) regardless of the range. On a d- 
dimensional (regular) network, Eq. Q still allows the analytic 
solution |* T , 17l3 ... 7 i d (i)| = J ni {2t) ■■■ J nd (2t)\, producing 
the behavior Pq(4) ~ t d . To characterize diffusion, we define 
the diffusion time r associated with the participation ratio by 
the condition 



P Q (t=r)=cN 



(8) 



with a constant c between zero and unity. For a d-dimensional 
network, we thus have the scaling behavior 



N 1/d . 



(9) 



For the WS network with shortcuts, the participation ratio in 
Eq. Q is computed numerically at given time t and the diffu- 
sion time t is measured from the condition in Eq. 0, where 
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FIG. 3: Diffusion time r versus the network size TV in the log-log 
scale. The case p = 0, corresponding to the local regular network, 
exhibits power-law behavior r ~ TV, whereas downward curvatures 
are present at nonzero values of p. Inset: r versus N for p = 0.5 in 
the semi-log plot: The behavior r ~ log N is manifested. Numerical 
errors are about the sizes of symbols. 

we choose the numerical factor c = 0.25. Other choice for 
the value of c does not make any qualitative difference in the 
scaling behavior of t, only if it is not too close to zero or unity. 

In Fig- we display the diffusion time r depending on the 
size N at various values of the rewiring probability p. Here 
the participation ratio Pq (t) has been computed for 3000 dif- 
ferent network realizations with given parameters N, r, and 
p, over which averages have been taken. It is observed that 
the dependence on N changes crucially as p is increased: At 
p = 0, the linear behavior r ~ N is indeed observed in agree- 
ment with the analytical result while the inset manifests the 
dependence r ~ log N forp = 0.5. We have not systemati- 
cally investigated the scaling behavior at smaller values of p; 
this is difficult from the computational point of view since 
N needs to be increased much to avoid finite-size effects. 
Nevertheless it is very plausible to conclude that the behav- 
ior t ~ log N persists as far as p is nonzero. The observed 
change of the scaling behavior from r ~ N to log N is very 
interesting in comparison with the change of the characteris- 
tic path length from / ~ N to log N already observed^ This 
suggests that the small-world transition at p = 0, separating 
the large-world behavior (I ~ N) from the small-world be- 
havior (I ~ log N), is also the fast-world transition between 
the slow-world behavior (r ~ N) and the fast-world behavior 
(t ~ log N) from the dynamical point of view. Note also that, 
in view of Eq. (|9j for a d-dimensional system, such logarith- 
mic behavior apparently indicates that the effective dimension 
of the WS network is infinite, which is consistent with the ob- 
servation of the mean-field nature. 4,5 

In Fig. [4] the diffusion time r versus the rewiring proba- 
bility p is displayed for various network sizes. As the size 
N becomes larger, the region described well by the relation 
t ~ P ^ covers a broader range of p, suggesting that the 
power-law behavior r ~ p^ 1 is valid at any nonzero value 
of p in the thermodynamic limit. The monotonic decrease of 
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FIG. 4: Diffusion time r versus the rewiring probability p in the WS 
network of various sizes. The behavior r ~ is observed in a 
range of p near p = 1, which expands as the size N is increased. 



t with p is easily understood: The more shortcut end points 
exist, the faster the diffusion is, as discussed above. Such 
anomalous quantum diffusion, observed in this work for com- 
plex networks, has also been investigated in various quasiperi- 
odic quantum systemsi^ii as well as in quantum systems 
which have classically chaotic counterparts 

For comparison, we also consider briefly diffusion of a clas- 
sical particle on the WS network. The particle is put on a 
randomly chosen vertex of the WS network constructed as be- 
fore, then allowed to hop to one of its neighboring vertices, 
chosen randomly at each time. While the return probability in 
such classical diffusion was examined^ it is appropriate here 
to consider, by analogy with the quantum mechanical partici- 
pation ratio Pq (t), the number of vertices visited by the parti- 
cle during time t, denoted by Pc(t). Note that this has the two 
limiting values: Pc(t=0) = 1 and Pc{t^oo) = N, which 
are the same as those of Pq. We thus call Pc the classical 
participation ratio for convenience. The diffusion time r for 
the classical diffusion is then determined by 

Pc(t=r) = O.liV, (10) 

where the use of numerical values other than 0. 1 again does 
not change the scaling behavior of r. In the absence of 
shortcuts (p = 0), the WS network reduces to the simple 
one-dimensional local regular network, where it is known 16 
Pc(t) ~ t 1 / 2 , and consequently we have the behavior r ~ 
N 2 . 

Figure [5] shows the diffusion time t versus the size N of 
the WS network for various values of the rewiring probabil- 
ity p. As in the quantum case, Pc (t) has been computed for 
3000 different network realizations and averages over them 
have been taken. As expected, for p = corresponding to the 
local one-dimensional lattice, the standard behavior t ~ N 2 
is observed. In sharp contrast, at any nonzero value of p the 
diffusion time displays the behavior r ~ N, demonstrating 
that classical diffusion on the WS network changes dramati- 
cally as soon as p takes a nonzero value. This behavior t ^ N 
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FIG. 5: Diffusion time r versus the network size N for classical 
diffusion. At p = the behavior r ~ A^ 2 is observed, whereas all 
curves at p ^ show the scaling behavior r ~ A" for sufficiently 
large AT. 
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FIG. 6: Classical diffusion time r versus the rewiring probability p 
for various network sizes. Unlike quantum diffusion, classical dif- 
fusion is not well described by the behavior of the diffusion time 



is also consistent with the infinite effective dimension of the 
WS network, in view of the known result that Pq (t) ~ t for 
dimension d > 3. 16 Comparing this observation for classi- 
cal diffusion with the previous one for quantum diffusion, one 
can draw in both cases the conclusion that the small-world 
transition at p = 0, associated with the change of the scaling 
behavior of the characteristic path length I, is accompanied 
by the fast-world transition in the dynamic view point, where 
the exponent a in the relaxation behavior r ~ A^" changes to 
a — 1. On the other hand, Fig.|5J where t versus p for clas- 
sical diffusion is plotted, shows that the behavior t ~ p^ 1 
observed for quantum diffusion (see Fig. |4j does not appear 
in classical diffusion. 

In summary, we have investigated both quantum diffusion 
and classical diffusion on Watts-Strogatz small-world net- 
works. To describe these, we have introduced the participa- 

TABLE I: Quantum and classical diffusion: the diffusion time r ver- 
sus the network size N. For comparison, the scaling behavior of the 
characteristic path length I is also presented. 

p = p / 

quantum r ~ N r ~ log N 

classical r ~ N 2 r ~ N 

characteristic path length Z ~ JV I ~ log N 



tion ratio and determine the diffusion time. The obtained scal- 
ing behaviors of the diffusion time, as summarized in TableU 
show dramatic changes as one increases the rewiring proba- 
bility from zero to a nonzero value. This is reminiscent of the 
small-world transition in the scaling behavior of the character- 
istic path length I from the small-world regime I ~ A^ to the 
large-world regime / ~ log N, and suggests that the small- 
world transition can also be termed the fast-world transition. 
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